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Abstract: 

The  average  refractive  bending  of  radio-frequency  rays  has 
been  calculated  from  actual  radio-frequency  refractive  index  dis¬ 
tributions  with  height  as  computed  from  mean  radiosonde  data.  The 
maximum  possible  bending  occurs  for  rays  passing  entirely  through 
the  atmosphere  and  arriving  or  departing  tangentially  at  the  earth' s 
surface.  The  range  of  this  total  angular  ray  bending  extends  from 
11  milliradians  (mr)  (0.63°)  at  Fairbanks,  Alaska,  in  April,  to 
18  mr  (1.01°)  at  San  Juan,  Puerto  Rico,  in  July,  and  is  about  14  mr 
(0.80°)  at  Washington,  D.C.,  in  October.  About  90$  of  the  ray  bend¬ 
ing  occurs  in  the  lowest  10  km  of  the  atmosphere.  These  results 
are  compared  with  the  ray  bending  computed  from  the  4/3  effective 
earth' s  radius  approximation. 
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*  This  paper  is  the  first  part  of  a  standard  refraction  study,  the 
goal  of  which  is  the  calculation  of  field  intensity  as  affected 
by  average  refraction  conditions  when  one  or  both  of  the  ter¬ 
minals  is  very  high,  such  as  a  high-flying  aircraft  or  rocket. 

By  the  "lower  atmosphere"  is  arbitrarily  meant  that  portion  of 
the  atmosphere  up  to  about  18  km.  However,  an  estimate  is  made 
in  this  study  of  the  contribution  to  the  refraction  of  the  at¬ 
mosphere  above  this  elevation,  excluding  any  ionospheric  effects. 
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1.  Introduction 

The  problem  of  radio-ray  tracing  ia  closely  connected  with  the 
problem  of  astronomical  refraction.  The  principal  difference  be¬ 
tween  optical-ray  and  radio-ray  refraction  is  due  to  the  effect  of 
the  large  contribution  of  water  vapor  to  the  refractive  index  of 
the  atmosphere  at  radio  frequencies.  The  history  of  atmospheric 
optical-ray  tracing  is  a  long  one  and  is  mostly  devoted  to  the 
computation  of  astronomical  refraction  and  somewhat  to  geodetic 
surveying.  Garfinkel'^  reviews  the  past  history  of  astronomical 
refraction  calculations,  and  he  presents  a  method  of  obtaining 
total  optical  refraction  for  a  ray  passing  entirely  through  the 
earth's  atmosphere  for  angles  of  elevation  less  than  10°.  In  all 
the  astronomical  refraction  Investigations,  the  air  density-height 
distribution  function  is  approximated  by  plausible  physical  models. 
Gar fink el  fits  mean  sounding  data  to  his  model,  which  consists  of 
a  temperature  lapse  up  to  the  stratosphere,  and  then  an  isothermal 
stratosphere.  He  shows  that  his  results  agree  quite  well,  within 
the  error  of  measurement,  with  the  sesai-emplrical  tables  of  Rad&u 
and  Pulkowa. 

The  tabulated  sounding  data  which  he  used  were  taken  from 
Humphrey f  ^  and  were  representative  of  four  European  cities  be¬ 
tween  the  years  1900  and  1912.  Since  the  date  of  Garfinkel's  paper 
(February  1944),  the  United  States  Weather  Bureau  has  published 
(May  1945)  averages  of  radiosonde  data  taken  for  the  most  part 
from  July  1939  to  December  1943  for  the  North  American  network  of 
stations  from  Alaska  to  the  Caribbean  Sea. '*)  These  data  are  used 
in  the  present  computations. 


(1)  "An  Investigation  in  the  Theory  of  Astronomical  Refraction,"  by 
B.  Garfinkel,  Astronomical  Journal,  Vol.  50,  No.  8  (2/44)* 

(2)  "Physics  of  the  Air,"  by  W.  J.  Humphreys,  McGraw-Hill  (1929) 
pp.  54^*56. 

(3)  "Upper  Air  Average  Values  of  Temperature,  Pressure,  and  Relative 
Humidity  over  the  United  States  and  Alaska,"  by  R.  Ratner, 
Climate  and  Crop  Weather  Div.,  U.S.W.B.,  Washington,  D.C. 

(5/45). 
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2.  Atmospheric  Refraction*  of  Radio  Wares. 

The  object  of  atmospheric  radio  wave  refraction  studies  usually 
is  to  find  a  convenient  method  of  correcting  the  radio  wave  field- 
intensity  equations  for  this  effect.  This  final  step  has  not  been 
taken  in  this  paper,  but  a  method  of  obtaining  the  true  refraction 
is  presented  in  Section  4,  page  8.  A  survey  of  existing  methods 
is  made  herein  in  terms  of  refraction  calculations,  instead  of 
field-intensity  applications  for  which  they  were  originally  de¬ 
signed,  and  it  is  shown  how  widely  divergent  the  results  of  the 
different  methods  are  from  the  true  refraction.  We  show  that  the 
problem  is  one  of  correcting  for  refraction  at  different  heights 
in  the  atmosphere  and  for  different  seasons  at  different  geograph¬ 
ical  locations.  The  most  widely  used  method  of  approximation  is  ob¬ 
tained  by  using  an  effective  earth's  radius  equal  to  4/3  the  actual 
radius.  We  furnish  data  to  examine  the  adequacy  of  the  4/3  earth's 
radius  approximation. 

In  the  present  paper,  refraction  has  been  computed  from  the 
actual  refractive  index  structure  of  the  atmosphere  as  obtained 
from  average  meteorological  data  at  the  location  in  question.  The 
assumption  is  made,  as  is  usually  done,  that  the  atmosphere  is 
horizontally  homogeneous  over  that  location.  Calculations  have 
bean  made  for  different  refractive  index  structure®  over  different 
locations.  The  results  have  been  plotted  aa  the  solid  curves  A, 

B,  C  in  Figure  1,  page  &•  Curves  A,  B,  and  C  represent  the 
actual  variation  of  refraction  with  height,  computed  by  methods 
presented  in  this  paper,  for  Fairbanks,  Alaska,  in  April,  for 
Washington,  D.C.,  in  October,  and  for  San  Juan,  Puerto  Rico,  in 
July.  Broken  curves  A* ,  B' ,  and  C1  represent  the  approximations 
obtained  by  a  method  due  to  Fearcey  (see  below)  using  a  4/3  earth's 
radius  factor  up  to  the  height  at  which  the  refractive  index  becomes 
unity  and  considering  no  refraction  beyond  this  point. 

The  earliest  successful  attack  on  the  problem  of  correcting 
radio  propagation  analysis  for  atmospheric,  refraction  was  intro¬ 
duced  by  Sohelleng,  Burrows,  and  Ferrell when  they  introduced 


*  Refraction  ia  defined  here  as  the  total  bending  of  a  ray  over 
a  path.  The  curvature  at  a  point  on  a  curve  ia  defined  as  the 
rate  of  change  of  the  tangent  direction  per  unit  length  of 
curve  at  that  point,  and  ia  measured  in  terms  of  change  of 
angle  per  unit  length  of  path.  For  a  circle,  the  curvature 
at  any  point  is  constant  and  ia  equal  to  the  reciprocal  of 
the  radius.  When  the  curvature  ia  known  along  a  ray,  the 
refraction  may  be  computed  as  the  integral  of  the  curvature 
over  the  path. 

(4) "Ultra  Short  Wave  Propagation,*  by  J.  C.  Sohelleng,  C.  R.  Bur¬ 
rows,  and  E.  B.  Farrell,  Proe.  I.R.E.,  Vol.  21,  p.  427  (3/33). 
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the  concept  of  effective  earth’s  radius.  They  found  that  a  factor  of 
4/3  was  satisfactory  for  annual  average  conditions,  up  to  0.5  km  as 
deduced  from  Humphreys’ '2)  weather  data.  However,  the  correction 
thus  introduced  is  only  strictly  correct  for  an  atmosphere  with 
linear  refractive  index  gradient  of  39  x  10“6  units  per  km,  and 
thus  cannot  possibly  hold  beyond  the  height  (approximately  10  km) 
at  which  the  refractive  index  has  reached  the  value  of  unity.  In 
addition,  the  atmosphere  has  a  non-linear  refractive  index  gradient, 
and  the  true  correction  will  vary  with  elevation  and  distance. 

Eckart  and  Plendl^  showed  how  to  correct  for  atmospheric 
refraction  for  nearly  horizontal  rays  by  fitting  the  actual  re¬ 
fractive  index  structure  by  means  of  two  parabolas,  one  from  the 
surface  to  5  km,  and  the  other  from  5  to  10  km.  Above  that  eleva¬ 
tion  they  considered  the  refraction  negligible.  The  method  fur¬ 
nished  satisfactory  tables  of  corrections  for  a  mean  atmosphere 
to  10  km.  To  extend  the  corrections  to  greater  heights  or  to  dif¬ 
ferent  refractive  index  distributions  as  encountered  in  different 
geographical  locations  would  require  additional  parabolic  fits  of 
data  and  extensive  calculations. 

(6  7) 

Stickland v  * '  attempted  to  correct  radio  field-intensity 
calculations  for  the  variation  of  refraction  with  height.  An 
empirical  relation  was  found  for  the  average  atmospheric  refrac¬ 
tive  index  distribution  in  terms  of  the  ratio  m  of  the  ray  radius 
of  curvature  at  any  height  h  in  the  atmosphere  for  nearly  hori¬ 
zontal  rays, 


m  *  2.5  +  1.0  h  (1) 

h  being  given  in  km.  Meteorological  conditions  in  the  lowest 
kilometer  are  admittedly  quite  variable,  but  the  formula  was  con¬ 
sidered  valid  to  10  km,  above  which  height  the  refraction  was 
assumed  to  have  ceased.  It  was  then  suggested  that  for  computa¬ 
tional  purposes,  an  average  ratio  m  be  used  between  any  two  points 
close  together  at  heights  h^  and  hg,  5  -  2.5  +  0.5(h^  +  l^) .  When 
one  point  is  near  the  ground,  h^  5  0,  and  then  we  have: 

(5)  "Die  Uberwindung  der  Erdkrummung  bei  Ultrakurzwellen  durch  die 

Strahlenbrechung  in  der  Atmosphere, "  by  G.  Eckart  and  H.  Plendl, 
Hochfreq.  md  Elektroakustik,  Vol.  52,  No.  2  (8/38)  pp. 44-58. 

(6)  "The  Effect  of  Atmospheric  Refraction  on  the  Propagation  of 
Radio  Waves,"  by  A.  C.  Stickland,  R.R.B./S.  10  (3/43). 

(7)  "The  Dielectric  Constant  of  Water  Vapour  and  its  Effects  on  the 
Propagation  of  Very  Short  Waves,"  ty  A.  C.  Stickland,  R.R.B./S. 2 
(5/42). 
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m  *  2*5  +  0.5  h  (2) 

If  we  use  this  average  m  for  calculating  the  actual  refraction 
of  a  tangential  ray,  we  get  incorrect  results.  In  order  to  demon¬ 
strate  this,  we  will  now  proceed  to  work  out  the  total  refraction 
in  terms  of  Stickland' s  average  m  concept,  and  show  that  these 
values  disagree  with  the  true  values  of  refraction  as  computed  by 
the  method  outlined  later  on  in  this  report.  The  Stickland  approxi¬ 
mation  calls  for  the  use  of  a  single  average  value  of  m  between  the 
ground  and  the  height  h  and  implies  that  we  may  use  for  this  layer 
that  value  of  effective  earth'  a  radius  factor  lc,  related  to  m  in 
equation  2  by* 

(earth's  curvature)  -  (ray  curvature)  *  ("effective"  earth's  curvature) 


1 

a 


1_ 

ns 


or.  E  =  (3) 

Now  a  general  expression  for  the  refraction  T  for  a  tangential  ray 
through  a  given  layer  of  depth  h,  with  an  effective  earth' s  radius 
factor  k  corresponding  to  some  linear  gradient  of  refractive  index  is: 

(See  Appendix  4»  page  20) 

where:  T  is  the  refraction  in  radians. 

a  is  the  earth's  radius,  expressed  in  the  same  units  as  h. 
h  is  the  depth  of  the  layer,  measured  from  the  ground, 
k  is  the  effective  earth's  radius  factor. 

The  bending  of  a  radio  ray  between  the  ground  and  height  h,  ob¬ 
tained  from  the  last  equation  by  using  the  values  of  k  obtained  from 
equations  (2)  and  (3)  has  been  plotted  in  Figure  2,  page  24»  Curve 
A.  In  this  same  figure  we  have  also  plotted:  (l)  Curve  B,  the  actual 
refraction  for  a  Stickland  atmosphere  defined  by  equation  (l)  and  com¬ 
puted  by  the  method  presented  in  this  report,  (2)  the  mean  refraction 
for  Washington,  D.C.,  in  October,  and  (3)  the  4/3  earth's  radius  ap¬ 
proximation.  It  may  be  seen  that  the  Stickland  averaging  process 
gives  too  low  values  of  refraction,  and  that  the  simple  4/3  earth' s 
radius  approximation  is  at  least  as  satisfactory  as  Stickland' s 
refraction  correction.  The  reason  for  the  failure  of  Stickland' s 
formula  in  this  case  is  that  the  method  of  averaging  m  has  the  ef¬ 
fect  of  weighting  all  height  intervals  equally,  whereas  the  concept 
of  the  refraction  as  the  integrated  curvature  along  the  ray  path 
allows  for  the  condition  that  almost  horizontal  rays  have  much 
longer  paths  in  the  lower  layers  where  the  refraction  is  greatest. 
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(8) 

Pearcey v  '  has  given  a  simple  arithmetical  formula  (equation  4 
below)  for  correcting  for  the  refraction  of  solar  radio  rays  arriving 
at  low  angles  0O,  a  procedure  necessary  in  solar  noise  studies,  by 
using  the  surface  value  of  the  refractive  index  and  then  assuming  a 
linear  refractive  index  gradient  of  39  x  1CP^  units/km  up  to  that 
height  at  which  the  refractive  index  becomes  unity.  Above  this 
height  the  refraction  is  considered  to  have  ceased.  This  is  equiva¬ 
lent  to  using  a  4/ 3  earth's  radius  factor  but,  in  addition,  the  final 
height  is  specified  at  which  the  correction  is  applicable.  This  model 
has  the  advantage  of  being  simple  and  the  validity  of  the  approxima¬ 
tion  for  total  bending  may  be  examined  in  the  dashed  curvesof  Figure 
1,  page  23.  Pearcey' s  formula  for  low  ground  angles  of  arrival  or 
departure,  0  £  80  =  10°  is  given  by: 


T 


Anp 


An.)  +  0 
o  o 


in  radians 


(4) 


where:  T  is  the  total  bending  in  radians. 

AnQ  is  (np-1)  where  rip  is  the  value  of  the  refractive  index 
at  the  earth's  surface. 

h2  is  the  elevation  at  which  An  becomes  sero  under  a 

standard  refractive  index  gradient  of  39  x  units/km. 

a  is  the  earth' s  radius,  expressed  in  the  same  units  as  h^. 


From  Figures  1  and  2,  it  may  be  seen  that  the  refraction  approxi¬ 
mations  in  use  at  the  present  time  are  not  quite  adequate  for  all 
seasons,  geographical  locations  or  at  very  high  altitudes.  However, 
the  effect  on  field-intensity  calculations  of  the  deviations  of  the 
refraction  approximations  from  the  true  refraction  has  not  been  de¬ 
termined  in  this  study.  It  is  contemplated  to  study  these  effects 
and  to  present  the  results  in  a  subsequent  paper.  It  is  also  in¬ 
tended  to  use  the  receivers  of  existing  radar  installations  to  get 
an  experimental  measure  of  the  total  atmospheric  refraction  of 
radio-frequency  noise  from  the  sun  at  one  or  two  frequencies  be¬ 
tween  100  Mc/s  and  10,000  Mc/s  for  comparison  with  the  refraction 
computed  from  nearly  simultaneous  meteorological  soundings. 


(8)  "Solar  Radiation  at  Radio  Frequencies  and  its  Relation  to 

Sunspots,"  by  J.  L.  Pawsey,  R.  Payne  Scott,  L.  L.  McCready, 

R -P.R.  24  (  6/46)  Commonwealth  of  Australia,  Council  for 
Scientific  and  Industrial  Research,  Division  of  Radiophysics. 
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3.  Meteorological  Data  and  Refractive  Index  Computations. 

The  formula  for  the  radio-frequency*  refractive  index  of  moist 
•ir  1.(9): 

(n-1)  *  106  =  22E  (1  +  4800  2-  )  (5) 

T  pT 

Equation  (5)  may  be  written  in  the  form: 

Optical  and  Radio  Explicit  Water  Vapor  Term 

Frequency  Term  Required  only  at  Radio  Frequencies 

(n-1)  x  106  =  79p/T  +  3.79  x  104  e/T2  (5a) 

where:  n  is  the  refractive  index. 

p  is  the  total  pressure  in  millibars  (mb)  and  is  equal  to 

Pd  +  •• 

is  the  partial  pressure  of  dry  air. 
e  is  the  water  vapor  pressure. 

T  is  the  temperature  in  °K,  and  is  equal  to  273  +  t°  C. 

This  is  a  semi-empirical  formula  and  was  obtained  as  a  result  of 
radio  wave  propagation  studies  and  laboratory  measurements  of  the  re¬ 
fractive  index  of  water  vapor  and  dry  air.  The  adequacy  of  equation 
(5)  as  a  basis  for  refraction  calculations  is  discussed  in  Appendix  1, 
page  12.  Charts  based  on  this  formula  have  been  constructed  for  ob¬ 
taining  the  refractive  index  of  moist  air  at  radio  frequencies  from 
meteorological  data  to  heights  of  25  km.  In  this  report,  representa¬ 
tive  radiosonde  data  published  by  the  U.S.  Weather  Bureau^ 3)  have  been 
used  for  the  calculation  of  the  refractive  index  distribution  for 
three  North  American  stations  separated  widely  in  geographical  loca¬ 
tion  for  one  month  of  each  season.  The  computations  were  made  by 

*  The  formula  is  applicable  for  frequencies  at  least  up  to  30,000  Mc/s. 

Although  some  variation  of  refraction  might  be  expected  to  exist 
near  various  absorption  lines  associated  with  various  elements  in 
the  atmosphere.  Van  Vleck'^'  has  pointed  out  that  the  observed 
magnitude  of  the  absorption  coefficient  of  oxygen  and  water  vapor  in 
.the  atmosphere  whose  absorption  lines  lie  below  30,000  Mc/s  would 
not  be  expected  to  be  accompanied  by  appreciable  changes  in  re¬ 
fraction.  Further  investigation  may  show  that  formula  5  is  ap¬ 
plicable  at  frequencies  even  considerably  in  excess  of  30,000  Me/e. 

(9) "Nomograms  for  Computation  of  Modified  Index  of  Refraction,"  by 
R.  A.  Burgoyne,  M.I.T.  RL  #551  (4/45). 

(10)  "The  Relation  between  Absorption  and  the  Frequency  Dependence  of 
Refraction,"  by  J.  H.  Van  Vleck,  M.I.T.  RL  #735  (5/45) 
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drawing  graphs  respectively  for  the  contribution  due  to  the  total 
pressure  term,  79  p/T,  important  both  at  optical  and  radio  frequencies, 
and  that  due  to  the  explicit  water  vapor  pressure  term,  3*79  x  10^  e/?  , 
important  only  at  radio  frequencies.  Figures  3-6,  pages  25-28 
give  the  contribution  of  the  total  pressure  term  to  the  refraotive 
index  computation.  Figures  7  and  8,  pages  29-30  ,  give,  in  terns  of 

relative  huaidity(RH)and  temperature,  the  contribution  due  to  the 
explicit  water  vapor  term  (Relative  HumidityJ’eoperature  Term).  The 
relation  between  e  and  RH  is  e  =  eg  (RH)/l00,  ea  being  the  saturation 
water  vapor  pressure  which  is  only  a  function  of  temperature .  The 
explicit  water  vapor  term  then  becomes  3.79  x  102  e®  (RH)  A2. 

Data  on  geographical  and  seasonal  distributions  of  refractive 
index  have  been  examined,  and  it  is  found  that  the  extremes  are 
represented  by  Fairbanks,  Alaska,  in  April  and  by  San  Juan,  Puerto 
Rico,  in  July.  The  data  for  Washington,  D.C.  in  October  were  found 
to  be  approximately  a  mean  between  these  two  extremes. 

Ordinarily  soundings  go  up  only  to  15  to  18  km.  However,  at 
these  elevations  the  atmosphere  is  approximately  isothermal  (strato¬ 
sphere),  and  the  refractive  index  distribution  may  be  easily  extended 
to  great  heights, at  least  to  30  k»,  which  is  as  far  as  the  isothermal 
condition  prevails^}.  (See  Appendix  2,  page  16). 

There  are  limitations  on  the  accuracy  of  both  the  refractive  index 
formula,  equation  (5a) ,  and  the  meteorological  observations  in  an  in¬ 
dividual  sounding.  These  limitations  are  discussed  in  Appendix  1, 
page  12,  where  there  is  shown  to  be  an  uncertainty  in  the  refractive 
index  calculations  of  the  order  of  7*£,  aboui  half  being  due  to  the 
uncertainty  in  the  formula  coefficients,  and  half  due  to  the  errors 
of  meteorological  measurement,  the  humidity  measurements  being  the 
least  reliable. 

4.  Computation  of  Refraction. 

A  method  is  now  presented  for  computing  the  refraction  of  radio- 
frequency  rays  through  an  atmosphere  of  known  refractive  index  dis¬ 
tribution,  under  the  usual  assumption  that  surfaces  of  equal  Index 
of  refraction  are  spherical  and, concentric  with  the  earth  (horisontal 
homogeneity) .  7h®  refraction  T*  ^  ^  of  a  s'ay,  as  illustrated  in  Figure 

9,  page  31,  between  two  heights  n{  and  %  through  an  atmosphere  of 
known  refractive  index  distribution  and  making  an  angle  ©  with  the 
spherical  refracting  surface  at  the  height  h  is  derived  in  Appendix 
3,  page  17,  ass 

T1j2  s  (cot  ^,2  ^  “  Aaj)  in  radian#  (6) 

Til)  "Tentative  Tables  for  the  Properties  of  the  Upper  Atmosphere".  - 

C.  N.  Warfield,  N.A.C.A.  Technical  Note  1200  (l/47). 
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where,  An  is  (n-l)  and  (cot  ©)«  2  the  mean  -value  ©f  the  cotangent 
function  over  the  interval  hj_,  ii.  Our  problem  is  solved  when  we  ®r® 
able  to  obtain  a  proper  value  or  (cot  ©)j  the  explicating  factors 
being  that  ©  is  a  function  of  both  An  ®M  h  as  determined  from 
Snell f  s  law  for  concentric  spherical  refracting  surface® ,  and  further 
that  An  is  only  known  empirically  as  a  function  of  h,  being  given  by 
a  table  of  measured  value®  as  obtained  from  radiosonde  data®  Thu® 
according  to  Snell8 e  law,  the  angles  ©  for  a  ray  path  through  con¬ 
centric  spherical  refracting  surfaces  are  determined  from  the 
formula, (l2) 

nr  cos  ©  «  cos  ©0»  constant  (7) 

where*  r  s  a  +  h,  is  the  radius  of  the  spherical  refracting  surface, 
at  the  elevation  h,  and  a  is  the  earth8®  radius,  and  where  the  sub¬ 
script  zero  refers  to  the  ground  level* 

At  small  elevation  angles,  @0  £  10°,  equation  (?)  transform®, 
neglecting  quantities  of  second  and  higher  orders,  t©vl2) $ 

*2  ~  8°  «  (  An  +  J)  -  <  An„  +  j*)  (8) 

and  we  then  have,  solving  for  ©, 

©  s  +  2{  An  +  |)  -  2(  An^  4  j*  )  (9) 

where  h^  is  the  station  elevation  above  mean  sea  level*  This  enables 
us  to  calculate  ©  at  any  elevation,  provided  that  the  calculated 
value  of  ©  is  not  appreciably  greater  than  10®,  since  we  are  given 
for  each  station  a  table  of  values  of  An  against  h*  The  An  data 
come  from  the  refractive  index  distribution  with  height,  which  are 
obtained  from  the  observed  meteorological  data,  as  discussed  in 
Section  3.  It  is  necessary  to  take  into  account  several  height 
intervals,  over  each  of  which  (cot  §7  may  b®  conveniently  deter® 
mined  in  the  application  of  equation  (6) *  It  is  shown  in  Appendix 
4,  page  19,  how  the  number  and  size  of  these  height  intervals  i© 
determined.  Only  four  height  intervale  are  required  from  the  ground 
to  18  km,  namely,  (l)  surfacA  to  0*5  kmj  (2)  0.5  to  2.5  km;  (3)  2.5 
to  6*0  km;  (4)  6  to  18  km.  Above  18  km  the  contribution  to  the 
total  refraction  is  small  and  may  be  shown  to  be  given  approximately 
by  T  *  Anig  cot  %g  where  the  subscripts  refer  to  the  values  of 
the  quantities  at  18  km.  Over  these  intervals  at  these  low  angles, 
it  may  be  shown  that  (cot  0)  s  1/6*,  where  6*  is  the  arithmetical 

(l2)  "Tropospheric  Propagation  and  Radio-Meteorology,"  WPG  5, 

Reprint  CRPL-T3  (10/46). 
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mean  of  ©  on  the  interval.  Thus,  for  the  interval  (hn,  h?),  we  have 
■  (©i  +  ©2)/2.  When  we  substitute  this  approximation  in  equation 
(6),  with  Qjg  in  radians,  we  obtain  simply, 

^  Ani  -  An*>  ^  < 

T  ~  —  ' ' "  — -  in  radians  (0  *  0Q  =  10°)  (10) 

t 

where  the  angles  6  are  determined  from  equation  (9)  above,  and  the 
An  are  obtained  from  meteorological  data  as  described  in  Section  3* 

Whenever  the  angle  of  arrival  or  departure  ©0  is  greater  them 
10°,  the  total  refraction  through  the  entire  atmosphere  is  very  small 
and  can  be  obtained  within  an  error  of  !&  from  the  simple  formula  de¬ 
rived  from  equation  (6)  in  Appendix  4,  page  19, 

T  s  An0  cot  ©0  in  radians  (10°  <  0Q  =  90°)  (11) 

the  approximation  implied  being  that  cot  0  =  cot  0O. 

The  mean  ray  bending  for  radio-frequency  rays  over  the  range  of 
angles  of  arrival  or  departure  has  been  obtained  for  San  Juan,  Puerto 
Rico  (July)  and  Fairbanks,  Alaska  (April),  as  examples  of  extremes} 
and  for  Washington,  D.C.  (October)  as  a  mean  case.  These  data  have 
been  plotted  in'  Figures  10  and  10a,  pages  32-33  •  In  these  same 
figures,  we  have  plotted  our  computed  values  of  optical  refraction 
for  Washington,  D.C.,  in  October  and  compared  them  with  the  mean 
observed  astronomical  refraction  data (13)  as  adjusted  to  these  same 
Surface  refraction  conditions.  Optical  refraction  for  Washington,  D.C. 
in  October  was  computed  by  the  methods  of  this  section,  using  the 
optical  refractive  index  distribution  obtained  by  neglecting  the 
explicit  water  vapor  term  in  the  radio-frequency  refractive  index 
computation.  Good  agreement  is  obtained  which  provides  a  check  of 
the  computational  procedure. 

/ 

It  is  now  possible  to  determine  the  ray  path  as  a  function  of 
the  distance  d  along  the  earth's  surface  as  illustrated  in  Figure  9, 
page  31.  ©  has  been  determined  from  equations  7  or  9,  and  2^*® 

been  determined  from  equations  10  or  11  and  it  may  be  shown  from  geo¬ 
metric  considerations  in  Figure  9  that: 

T  =  '  +  (®i  -  ©2)  in  radians  (12) 

and  we  then  have:  d^g  1  8  (  ^  1  2  +  ‘  ®1  ^  (13) 

where:  2*  are  expressed  in  radians,  and  d  is  expressed 

in  the  same’unit  as  a,  the  earth's  radius. 

(l3)  "Astronomy, "  Vol.  1—Solar  System  by  H.  N.  Russell,  R.  S.  Dugan 

and  J.  Q.  Stewart,  Ginn,  1945. 
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5.  Conclusion 

It  has  bee®  found  that  ifa®  actual  air  refraction  in  the  lower 
atmosphere  deviates  considerably  froa  that  which  has  been  assumed 
when  we  use  an  effective  radius  of  the  earth  equal  to  4/3  of  its 
actual  value  as  an  allowance  for  this  refraction,  Thu@?  although 
the  4/3  earth*  s  radius  aasisaptioa  provides  a  reasonably  g@©d  ap¬ 
proximation  in  the  average  eas@p  these  deviations  £r©a  the  true 
refraction  have  been  shown  to  fee  fairly  large  both,  as  a  function 
of  altitude  and  ©f  geographical  location,  Th@  practical  importance 
of  these  deviations  in  terms  of  field-intensity  calculations  has 
not  as  yet  been  evaluated ,  but  they  may  turn  out  to  require  another 
type  of  refraction  approximation  which  allows  for  the  observed 
variations  with  respect  to  height  and  geographical  location. 
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Appendix  1.  Estimate  of  Uncertainty  In  Application  of 

Refractive  Index  Formula  for  Moist  Air. 

There  are  limitations  on  the  accuracy  of  both  the  refractive  index 
formula  and  the  meteorological  observations  in  an  individual  sounding. 

A.  Uncertainty  in  the  Refractive  Index  Formula. 

The  formula  for  the  refractive  index  of  moist  air,  equation  (5), 
page  7,  is  based  on  a  series  of  separate  laboratory  determinations  of 
the  dielectric  constants  of  dry  air  and  water  vapor ) .  The  contri¬ 
butions  of  each  are  then  combined  to  yield  equation  (5), 

(n-l)  x  106  s  22E  (i  +  4800  e/pT) 

T 

(a)  Dry  Air  Refractive  Index 

The  refractive  index,  n^,  of  dry  air  is  given  by: 

(nd  -  1)  x  106  -  79pd/T  (16) 

where  pd  is  the  partial  pressure  of  dry  air.  There  have  been  many 
measurements  taken  on  the  dielectric  constant  of  dry  air.  The  most 
recent  determination  has  been  made  by  Hector  and  Woernley(lA) .  From 
their  data  the  coefficient  in  equation  (16)  is  computed  to  be  76.5, 
whereas  the  average  value  obtained  by  all  workers  in  the  field  is 
77.9.  It  is  pointed  out  in  the  same  paper  that  the  internal  consis¬ 
tency  of  each  worker' s  data  has  been  better  than  the  agreement  from 
observer  to  observer.  The  entire  range  of  the  sixteen  determinations, 
the  first  of  which  was  made  by  Boltzmann  in  1874,  is  from  72.8  to 
80.8.  The  computed  probable  error  is  4^  of  the  mean.  It  should  be 
pointed  out  that  the  coefficient  in  equation  (16)  is  about  larger 
than  the  mean  value  for  all  workers  in  the  field.  Hence,  we  take  5^ 
as  a  reasonable  value  of  the  uncertainty  in  the  dry  air  contribution. 

(b)  Water  Vapor  Refractive  Index. 

The  refractive  index,  ny,  of  water  vapor  is  given  by: 

(nv  -  l)  x  106  =  22®  (i  +  48OO/T)  (17) 

According  to  reference  9,  the  best  determinations  over  a  range 
of  temperature  values  are  those  due  to  Sfinger  and  those  due  to 
Stranathan: 

f  %  (68  +  [  X  105  )  singer  (18) 

,  )  T  T 

(a*  -1)  x  106  =  S 

/  ft  (73  +  2m 22-Z-lQl  )  Stranathan  (19) 


(14)  "Dielectric  Constants  of  Eight  Gases,"  by  L.  G.  Hector  and  D.  L. 
Vtoemley,  Phys.  Rev.,  Vol.  69,  Nos.  3,  4,  p.101  (2/46) 
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The  value  of  the  water  vapor  refractive  index  actually  used  in 
obtaining  equation  (5)  is  given  by  equation  (17): 

(iiv  -  i)  X  106  =  Z2s  (1  +  4SQ0) 


instead  of  the  mean  value  of  equations  (18)  and  (19) ,  which  is: 

(uy  -  1)  x  106  =  22«  (o#89  +  ^252  ) 

T  T 

Even  assuming  a  temperature  of  50°C,  the  values  of  (n^  -  l)  x  10^ 
obtained  from  equation  (17)  do  not  differ  from  those  calculated  from 
equations  (18)  and  (19)  by  more  than  2.5^. 

We  consider  now  the  uncertainty  in  the  formula  of  equation  (5) 
for  (n  -  l)  *,10”  for  moist  air  due  to  the  respective  errors  in 
(n^j  -  l)  x  10”  for  dry  air  and  (nw  -  1)  x  10”  for  water  vapor.  The 
surface  values  of  these  quantities  for  a  typical  station,  like  Wash¬ 
ington,  D.C.,  in  October  is  (n-l)  x  10&  =  332  units,  which  is  composed 
of  (ng-l)  x  10^  =  281  units  and  (ny-l)  x  10®  =  51  units.  Thus  a  5 % 
error  in  the  dry  air  contribution,  and  a  2.556  error  in  the  water  vapor 
contribution,  lead  to  a  weighted  error  in  the  formula  for  the  refractive 
index  of  moist  air  of  455. 

B.  Uncertainty  in  Meteorological  Data. 

(a)  Drv  Air  Refractive  Index 

The  following  error  data  used  in  discussing  the  dry  air  refractive 
index  has  been  computed  from  data  in  U.S.  Weather  Bureau  memoranda 
dated  May  15,  1943  and  May  19,  1943,  and  is  based  on  Gregg' s  Standard 
Atmosphere  which  represents  average  U.S.  data  at  40°N  latitude. 

Table  1:  Probable  Qrror  in  Computed  Pressure  and  Temperature  at  Fixed 
Levels  for  a  Single  Sounding  for  U.S.  Weather  Stations. 


Fixed 

Altitude 

Pressure 

mb 

%  Error  in 
Computed 
Pressure 

Temp. 

°C 

%  Error  in 
Computed 
Tempera ture 

%  Error 

(®d-l)  X 

10,000  ft 

701.1 

to.i 

-  1.4 

±0.3 

±0.4 

20,000  ft 

471.8 

±0.28 

—20. 6 

±0.48 

±0.8 

10  km 

270.5 

±0.59 

-44.5 

±0.66 

±1.3 

13  km 

171.0 

±0.82 

-55 

±0.69 

±1.5 

16  km 

107.0 

±0.93 

-55 

±0.69 

±1.6 
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The  error  below  altitudes  of  20,000  f t  (  ~  6  km)  in  (n-l)  x  10^ 
for  moist  air,  due  to  the  error  in  the  dry  air  refractive  index,  must 
be  less  than  about  0.8^.  At  10  km  and  higher,  the  error  is  larger, 
but  there  is  a  negligible  water  vapor  contribution  at  these  altitudes 
and  so  the  tabulated  error  then  represents  the  total  error.  Belov 
about  10  km,  we  must  also  consider  the  effect  of  the  water  vapor 
refractive  index  contribution,  (ny-l)  x  10°. 


(b)  Water  Vapor  Refractive  Index 


For  our  purposes  the  relative  humidity  error  is  taken  as  ±10£  RH 
(  ±10  RH  units)  at  all  values  of  relative  humidity.  At  elevations 

from  10  km  on  up,  the  water  vapor  contribution  is  negligible  because 
of  the  negligible  amount  of  water  vapor  at  these  elevations.  The 
water  vapor  contribution,  equation  (17),  may  now  be  written: 


(V-l)  x  lo1 


,6  79es  (RH) 


(17) 


where  the  last  approximation  holds  at  ordinary  temperatures.  Now,  we 
see  in  Table  1  that  the  error  in  the  temperature  measurement  is  of  the 
order  of  1%  below  10  km.  An  examination  of  equation  (17)  will  show 
that  the  error  in  (ny-l)  x  10°  for  water  vapor,  due  to  the  tempera¬ 
ture  error,  is  of  the  order  of  2 t.  The  corresponding  error  in  e8, 
the  saturation  specific  humidity,  due  to  the  1^  error  in  T,  is  less 
than  1%.  The  error  due  to  the  uncertainty  in  the  relative  humidity 
measurement  has  been  calculated  for  a  typical  station,  taken  as . 
Washington,  D.C.,  in  October,  and  tabulated  in  Table  2,  page  15* 

The  total  percentage  error  in  (ny-l)  x  10°  for  water  vapor,  listed 
in  column  3,  was  obtained  by  adding  a  1*£  error  due  to  ea,  a  error 
due  to  the  effect  of  the  temperature,  and  the  equivalent  percentage 
error  due  to  an  absolute  error  of  ±10  RH  units  for  any  value  of  rela¬ 
tive  humidity.  At  the  surface,  the  error  data  are  omitted  sinoe 
humidity  measurements  are  made  with  a  sling  psychrometer  and  the 
accuracy  of  the  data  far  exceeds  that  obtainable  with  the  radiosonde 
humidity  element. 

We  thus  see  that  there  is  a  weighted  percentage  error  of  about 
3^  below  20,000  ft  (  £  6  km)  in  the  refractive  index  computation  of 
moist  air  due  to  errors  in  measurement  of  the  meteorological  elements. 
The  3^  includes  2^  for  the  water  vapor  contribution,  and  1^  for  the 
dry  air  contribution. 

Our  conclusion  then  is  that  the  total  uncertainty  in  the  refrac¬ 
tive  index  calculations  is  of  the  order  of  7^,  about  being  due  to 
the  uncertainty  in  the  coefficients  of  the  formula,  and  about  3%  due 
to  the  errors  of  measurement  of  the  meteorological  factors. 


Table  2:  Error  in  (a, -l)  x  10^  and  Corresponding  Error  in  (n-l)  x  10^  for  Moist  Air 
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Appendix  2.  Extrapolation  of  Refractive  Index  Data  to  30  km 


The  expression  for  the  pressure  In  an  Isothermal  atmosphere  Is 
given  by: 


-(h-ho)/R'Ts 

P  =  P0« 


(14) 


uheret  Tg  is  the  constant  temperature  of  the  stratosphere  in  °K, 
but  which  may  vary  with  season. 


p^  is  the  surface  pressure  in  mb  which  in  this  case  refers 
to  the  base  of  the  stratosphere,  but  the  actual  value  is 
unimportant  for  the  present  application  as  shown  below. 

R'  is  the  gas  constant  for  the  dry  stratosphere. 


h-ho  is  the  height  referred  to  the  base  of  the  stratosphere. 


Assuming  no  water  vapor  at  these  elevations,  the  refractive  index 
is  given  by: 


(n-l)  x  106  *  79p/T8  s 


#-(h-ho)/R‘T. 


and  hence,  log(n-l)  *  k^  -  l^h 


(15) 


where  and  kg  are  constants.  Thus  by  plotting  the  stratosphere  data 
that  we  have  for  the  location  concerned  (see  Figure  11,  page  34),  we 
should  get  a  straight  line  which  allows  us  to  extrapolate  refractive 
index  data  to  great  heights.  In  fact,  we  find  that  (n-l)  reduces  to 
4  x  ICH®  units  at  an  elevation  of  30  km  at  Washington,  D.C.  in  October, 
if  we  assume  the  stratosphere  to  be  isothermal  to  that  elevation.  'H' 
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Appendix  3i  Derivation  of  Formula  for  Refraction.  T  .  of 

a  Radio  Ray. 

Snell's  lav  for  spherical,  refracting  surfaces  is:  nr  cos  9  *  constant 
The  expression  for  curvature  at  a  point  on  the  ray  path,  to  be  derived 
below  from  equation  (7),  is: 

il.  .  i  «  -  Ida  cos  9  (20) 

4s  P  n  dh 


where:  P  is  the  radius  of  ray  curvature  at  a  point  on  the  ray. 

(— )  is  the  vertical  gradient  of  refractive  index  at  the  point  in 
question, 

/  dXx 

(  r*j”)  is  the  rate  of  change  of  curvature  at  any  point  along 


the  path. 


Differentiating  equation  (7)  (shown  above),  with  respect  to  height, 
h,  vemenbering  that  r  =  a  +  h,  we  have* 


1  4fi+  1  -  tan  9^*0 
n  dh  r  dh 

and 

r  008  e  “  8in  6  iS  *  *  n  I?  cos  © 

(21) 

From  the 
we  have: 

infinitesimal  section  of  the  ray  path  in  Figure  12,  page 

35» 

ds2  ■  dh^  +  40^ 

and 

r  *  cos  9 

ds 

and 

4h  »  sin  6 
ds 

and,  substituting  in  equation  (21),  for  cos  0  and  sin  9,  we  have: 

UMI  s  .  i  4a  cos  e 

ds  n  dh 

From  geometrical  considerations  for  the  quadrilateral  indicated  by 
dashed  lines  in  Figure  12, 


(22) 


( J|  +  8)  4.  (  ^  .  At  )  + 


and  hence,  At  *  A0  -  A0 


*0  +  [£•  -  (e 


+  A0) 


s  2'rr 


and, 


JX 

As 


a(0  -  e) 

As 
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Passing  to  the  limit  as  As  0,  ve  have,  substituting  in  aquation  (22), 

|X_  .  .  1  as  ooa  e  (20) 

4s  n  dfa 


Integrating  between  the  altitudes  h^  and  h^, 

=  1  ^  cos  e  da  ■  -f  2 

and  since  n  ■  1,  ^"1,2  ^  co^  ®  dn 


ig  -  ^ 


^  and  changing  the  variable  of  integration  to  An,  by 

means  of  the  relation  n  =  1  +  An 

An2 

wa  have  ^*1  2  =  ®  *(  An) 

^  An^ 

By  the  definition  of  the  mean  value  of  cot  6  over  the  interval 

(  Anj.,  Anj),  ^Ad, 


(cot  ®)i>2 


A®x 


cot  9  d(  An) 


(  An 2  -  An^) 
and,  therefore,  ve  have  finally, 

^1,2*  rCOt0)1>2  (Aai-ABj) 


(23) 


(6) 


Appendix  4.  Applications  of  Formula  for  Refraction 

of  a  Radio  Ray. 


To  determine  the  smallest  number  and,  therefore,  the  largest  sis® 

of  the  height  intervals  needed,  at  angles  10°,  consistent  with  the 
7t  error  arising  from  the  computation  of  An  as  shown  in  Appendix  1, 
page  12,  we  consider  the  following  facts: 


Under  the  condition  that  — 

dh 


An-  An-j 

”h  -"h. 


is  constant  in  the  in¬ 


terval  (h2,  h^)  we  have. 


Tcot~e)1  2  = 


L 


An. 


cot  8  d(  An) 


L 


cot  8  dh 


An2  -  An^ 


h  "  hl 


Now,  under  the  same  condition,  Snell’s  law  for  low  angles  of 
arrival  in  the  interval  becomes, 

9  “  ®1  „  ( An  +  J  )  -  (  Anx  +  ^  +  ^)(h  -  hj 

— 5  a  1  a  h  -  h,  a  - 


where 


h  -  h] 
ka 


(24) 


We 


ca  2  a  -5  >  P  being  the  radius  of  curvature  of  the  ray. 

2  r 

This  shows  8  to^be  a  linear  function  of  h  over  the  interval, 
therefore  plot  8^  against  h,  to  examine  how  wide  are  the  intervals 
over  which  this  linearity  is  indicated.  An  example  of  this  relation¬ 
ship  is  shown  in  Figure  13,  page  36, for  the  case  80  »  0  for  Washington, 
D.C.,  in  October.  This  empirical  relation  makes  it  possible  to  intro¬ 
duce  a  great  simplification  in  determining  the  value  of  (cot  8)  over 
such  an  interval.  For,  at  these  low  angles  cot  9  *  1/8,  and  it  may  be 
shown  that  cot  8  ~  l/@m,  where  is  the  arithmetical  mean  of  8  on  the 


interval,  e.g. ,  over  the  interval  (hi,  h2) ,  ©m 
Snell' 8  law  nrcos  ©  “  constant,  we  have: 


(81+  82)/2. 


From 


.  ,  d© 

cot  ©  =  ka 


Thus, 


_ h 


cot  0  dh 


cot  ©  - 


h2  -  hx 


ka 


92  “  gl 

h2  -  hx 


r  1 

"5 


and  therefore,  for  small  angles,  cot  8  s  ^  ^  ^ 

1  2 

When  this  value  of  cot  8  is  substituted  in  equation  6,  we  obtain,  simply, 

•■'i  -  ^ 


1,2  * 


in  radians 


0  =  @0  I  10° 


In  the  example  shown  in  Figure  13,  we  see  that  ©  can  be  considered 
arbitrarily  to  be  a  linear  function  of  h  over  four  separate  intervals  from 
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the  ground  to  18  km,  namely,  (l)  surface  to  0.5  km;  (2)  0.5  to  2.5  km;  (3) 
2.5  to  6.0  km;  and  (4)  6  to  18  km.  Above  18  km,  the  contribution  to  the 
total  refraction  is  small  and  may  be  shown  to  be  given  approximately  by 
T  »  cot  ,  where  the  subscripts  refer  to  the  values  of  the 

quantities  at  18  km.  In  general,  tlje  same  intervals  may  be  used  at  all 
locations,  since  the  shape  of  the  6^  vs  h  curve  does  not  change  appre¬ 
ciably  from  location  to  location. 

We  show  at  this  point  a  simple  formula  for  the  refraction  of  a  radio 
ray  tangential  to  the  earth' s  surface  in  an  atmospheric  layer  in  which 
an  effective  earth's  radius  factor  k  may  be  used.  Equation  (24)  shows 
that  in  this  case, 


In  addition  the  well-known  distance-height  formula  for  a  ray  tangent 
to  the  earth' s  surface  in  an  atmosphere  with  effective  earth’ s  radius 
factor  k  is: 

d  =  >/2ka  h 

Substituting  these  relations  in  equation  (12) ,  page  10,  we  have, 

tl,2*^*2+  <12> 

and  X » (  y/k  ) ,  ©x  being  0. 

At  large  angles,  0-  10°,  we  now  derive  the  following  expression 

for  the  total  refraction, 

T  s  Ano  cot  0O  in  radians  (ll) 

the  approximation  being  implied  that  (cot  §)  2  cot  ©0,  in  equation  (6) 

2  =  (cot  0)^2  ^  Anl  “  An2^  (6) 

We  use  the  subscript  o  for  the  ground  value,  the  subscript  18  at  h  *  18  km, 
owl  the  subscript  bo  for  h  >  a  Equation  (6)  may  then  be  transformed  to: 

X  T*  Teot  e)0>18  (  An.  -  A^g)  +  («>t  e)18(<D(  A"18  *  °>  <25> 

T  .  (cot  §T0(18  4n0  -  An18  [(cot  e^g  .  (cot  S)18jCC] 


or 


(25a) 
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For  average  refractive  index  distributions,  i.e.,  there  are  no 
height  intervals  where  the  downward  ray  curvature  exceeds  earth  cur¬ 
vature  and  neglecting  ionospheric  effects,  then  as  h  increases.  An 
approaches  zero  and  0  becomes  uniformly  larger  because  of  the  geometry 
of  the  ray  path.  This  means  that  cot  0  decreases  uniformly  with 
height,  and  hence,  the  average  value  (cot  ^  is  less  than 

(cot  0)o  Now,  Anig  is  l®38  than  An^lO  at  all  geographical 

locations  and  for  all  seasons. 


(cot  §y  . 

Also,  for  angles  0Q  a  10°  or  higher,  x -* — ~  - 

T55r5To,18 


Therefore,  at  most  the  second  term  in  equation  (23a)  is  5^  of  the 
first  term. 


In  addition,  from  consideration  of  data  at  San  Juan,  Puerto  Rico, 
in  July,  using  the  large-angle  fomula  (eq.(ll))  for  0  «  10°  which  is 
the  case  of  largest  refraction  considered  in  this  report,  the  largest 
error  we  make  by  taking  (cot  ©)Q  5  cot  0O  is  Thus,  at  0Q  *  10°, 

the  error  made  using  the  high-angle  formula  is  of  the  order  of  9%, 
with  the  error  decreasing  rapidly  as  Qq  increases  from  10°. 

Table  3  shows  a  sample  computation  using  the  ray  tracing  method 
of  this  report. 


22 


b3 

<j 

w 


„  gu 

o 


c 

O'  I  rn 1 

^  7  x 

•H 

(3 
<3 


<*> 

+ 


80  i— I 

« 

•H 


-ri  U 

r-  ®  s 


CM 

nO  *H 

© 


CM 

- — 

JECM 

«tn  nfjH 

O  II 


CM 


5E 


'© 


ON 


S3 

<a 


'o 


CM 


x 

as 

£■ 


CM 


On 

& 

to 


ON 

r-S 


CM 


flO 

CM 


CM 

CM 


•  a  s  s  s  r 


CV 


c*“\ 


o 

r-4 


8 

H 


ON 


ON 

r? 


vO 

• 

rH 

ON 


O' 

• 

W 

ON 


UN 

«r\ 


CM  O  CM 

*  £  £ 

H  UN 


UN 


£  $  $  S 

INI  f-  UN 


vO  O'  <M 

ON  «0  ON 

ON  ON  vO 


CM  © 
ON  H 
ON  ON 


& 

a 


s 


CM 

UN 


CM 


* 


£  $  3  R 

ON  O'  ^0 


o 

* 


UN 

CM 

O 


UN 


CM 


O 

vO 


to 


o 

00 


*3 

*r 

£ 

*8 

r-l 

■3 

o 

I 

8 


■a 


®  • 

■H  ® 

5 

13 

CO  o 


Tj  H 
*  H 

5* 

®  S3 
t)  *H 

x 

a  • 

S  ■? 

O  S3 1 


®  £ 

■P  2 

§£ 

rH  © 

i 

51? 

a  «h 

►»*d 


•  a 

£  S 

o 

a  ' 


5 

+ 

S3 

< 

a 

x 

a 

s« 

o 


REFRACTION,  T,  DEGREES 
0.4°  0.5°  0.6°  0,7°  0.8°  0.9 


133d  dO  SQNVSnOHl  V.1H9I3H 


23 


O' 


UJM*M*1H9I3H 


BENDING  BETWEEN  GROUND  AND  HEIGHT,  h.,  OF  A 
RADIO  RAY  TANGENTIAL  WITH  EARTH'S  SURFACE. 


REFRACTION,  T,  (DEGREES) 

0.2°  0.3°  0.4°  0.5°  0.6°  0.7°  0.8°  0.9°  1.0°  1.1°  1.2°  1.3° 
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Fig.  3  pressure -Temperature  term  for  the  refractive  index  computation 

(NUMBERS  ON  SLANT  LINES  ARE  VALUES  OF  79ft 
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TERM  FOR  THE  REFRACTIVE  INDEX 
SLANT  LINES  ARE  VALUES  OF  79^) 
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Fig.  5  PRESSURE-TEMPERATURE  TERM  FOR  THE  REFRACTIVE  INDEX  COMPUTATION 


(NUMBERS  ON  SLANT  LINES  ARE  VALUES  OF  79^) 
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Fig.  6  PRESSURE -TEMPERATURE  TERM  FOR  THE  REFRACTIVE  INDEX  COMPUTATION 

(NUMBERS  ON  SLANT  LINES  ARE  VALUES  OF  79^) 
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LEGEND 


=  EARTH'S  RADIUS 


= 


FATH  OF  RADIO  RAY,  IN  THIS 
CASE  LEAVING  THE  EARTH. 

ANGLE  OF  RADIO  RAY  WITH  THE 
EARTH'S  SURFACE  (ANGLE  OF 
ARRIVAL  OR  DEPARTURE ) 

ANGLES  OF  RADIO  RAY  WITH  THE 
SPHERICAL  REFRACTING  SURFACES 
AT  HEIGHTS  h,  AND  hg RESPECTIVELY. 

PROJECTED  DISTANCE  AUONG  EARTH'S 
SURFACE  OF  RAY  PATH  BETWEEN 
HEIGHTS  h,  AND  h2. 

REFRACTION  OR  ANGULAR  BENDING 
OF  RAY  BETWEEN  HEIGHTS  h,  AND  hg. 


Fig.  9  REFRACTION  T,  2  IN  THE  EARTH'S  ATMOSPHERE 
BETWEEN  HEIGHTS  h,  AND  h2 


ANGLE  OF  ARRIVAL,  0O,  (DEGREES) 

1°  2°  3°  4°  5®  6°  7°  8°  9°  10° 


:  -0.9® 


A-RADIO  RAY  REFRACTION,  SAN  JUAN,  PUERTO  RICO  (JULY) 

B-RADIO  RAY  REFRACTION,  WASHINGTON,  D.  C.  (OCTOBER) 

C-RADIO  RAY  REFRACTION,  FAIRBANKS,  ALASKA  ( APRIL) 

D-  ©  OPTICAL  REFRACTION,  WASHINGTON,  D  C.  (OCTOBER) 

COMPUTED  FROM  OPTICAL  REFRACTIVE  INDEX  OBTAINED 
FROM  AVERAGE  METEOROLOGICAL  DATA. 

A  MEAN  OBSERVED  ASTRONOMICAL  REFRACTION  [RUSSELL,  DUGAN,  STEWART  (1 3)] 
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Fig.  II  REFRACTIVE  INDEX  DISTRIBUTION  IN  THE  TROPOSPHERE  AND  STRATOSPHERE 

( WASH  I NGTON  ,0.0.,  OCTOBER ) 
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Fig.  12  INFINITESIMAL  SECTION 
OF  RAY  PATH 
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Fig.  13  <9*  VS.  HEIGHT  (<9  IS  THE  ANGLE  OF  ARRIVAL  AT  HEIGHT  h) 

WASHINGTON,  D.  C.,  OCTOBER,  <90  =  0  mr. 
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